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List  of  Acronyms  /  Abbreviations 


AM  Amplitude  Modulation 

AM-DSB-SC  Amplitude  Modulation,  Double  Side-band  Suppressed  Carrier 

ARMA  Auto-regressive  Moving  Average 

CBPF  Complex  Band-pass  Filter 

DAC  Digital-to- Analog  Converter 

DSP  Digital  Signal  Processing 

FIR  Finite  Impulse  Response 

FM  Frequency  Modulation 

I  In-phase 

HR  InBnite  Impulse  Response 

LPF  Lovy-pass  Filter 

NBIC  Narrow-band  Interference  Canceller 


Q 


Quadrature-phase 


1.  Introduction 

The  concept  of  the  Hilbert  transform  and  the  related  application  of  analytic  signal 
representation  are  well  known  and  studied  [1].  It  is  common  practice  in  modem  communications  and 
radar  signal  processing  systems,  particularly  those  that  are  implemented  digitally,  to  leverage  the 
benefits  of  analytic  representation.  Among  the  more  important  benefits  are  the  potential  increase  in 
spectral  resolution  for  a  given  number  of  signal  samples,  and  the  ability  to  calculate  instantaneous 
envelope  and  phase  measurements  from  the  resulting  complex-valued  signal.  What  is  less  well 
known  is  the  fact  that  there  are  alternative  methods  of  analytic  signal  generation  that  can  alleviate 
problems  associated  with  implementation  issues  such  as  real-time  operation  and  signal 
synchronization.  The  goal  of  this  report  is  to  bring  to  light  the  implications  of  recent  publication  in 
this  area,  and  to  augment  these  results  with  recent  AFRL  in-house  research  results.  In  particular, 
recent  published  research  introduces  filter  design  considerations  which  can  reduce  computational 
complexity  and  simplify  the  design  process  [2].  In  this  report,  this  method  is  generalized  for  a  wider 
range  of  filter  designs  and  applications. 

A  quick  review  of  the  Hilbert  transform  and  analytic  signal  representation  is  given  in  Section 
1.1.  For  more  detail,  the  reader  is  referred  to  [1].  Herein  we  are  concerned  with  digital  (i.e., 
numerical),  uniformly  sampled  discrete-time  sequences  resulting  from  the  proper  acquisition  of  a 
corresponding  analog  signal  of  interest  or  fi'om  the  computer  generation  of  a  signal  of  interest.  In 
practice,  one  can  certainly  approximate  Hilbert  transforms  and  the  analytic  signal  representation  with 
analog  waveforms.  This  method  is  described  in  Section  2.  However,  it  can  be  shown  that  even  with 
rather  stringent  specifications  on  in-phase  (I)  and  quadrature  (Q)  channel  balance  (e.g.,  0.1  dB 
amplitude  and  1  degree  phase  imbalance),  a  modest  amount  (40  dB)  of  image  distortion  suppression 
will  be  achieved  [3].  In  many  cases,  the  in-phase  and  quadrature  channels  are  then  digitized  in 
preparation  for  further  digital  processing.  In  contrast,  when  generating  the  analytic  signal 
computationally  after  digitization  of  the  signal  of  interest  [4,  5],  the  precision  of  the  required 
calculations  and  the  available  processor  resources  (speed,  memory,  etc.)  become  the  limiting  factors 
on  performance.  This  second  method  is  described  in  Section  3,  along  with  advancements  to  the 
method.  Because  of  the  recent  advances  in  digital  signal  processing  (DSP)  hardware,  this  latter 
approach  is  often  deemed  more  cost  effective,  particularly  when  image  suppression  greater  than  40  dB 
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is  required.  Section  3  addresses  often  overlooked  synchronization  concerns,  and  provides  a  novel 
interference  cancellation  application.  Section  4  provides  a  summary  and  conclusions. 

1.1  Review  of  Analytic  Signal  Representation 

A  convenient  method  of  explaining  the  concept  of  analytic  signal  representation  relies  on 
knowledge  of  some  basic  signals  and  systems  theory.  Although  the  term  “analytic”  has  roots  in  the 
theory  of  complex  variables,  it  will  suffice  here  to  say  that  this  term  refers  to  complex-valued 
continuous-time  /  continuous-amplitude  signals,  having  specific  continuity  properties  within  the  I  and 
Q  channel  signal  components.  In  this  sense,  it’s  somewhat  of  a  misnomer  to  refer  to  a  digitized  signal 
as  being  analytic,  given  the  discrete  nature  of  the  signal.  One  can,  however,  generate  good  analog 
approximations  to  a  true  analytic  signal  by  using  a  pair  of  digital-to-analog  converters  (DACs)  and 
low-pass  filters  using  the  corresponding  digital  “analytic”  signal.  Therefore,  for  the  remainder  of  this 
report  we  dismiss  distinctions  between  the  two,  allowing  the  term  “analytic  signal”  for  the  analog  and 
digital  cases.  Herein  we  are  ultimately  interested  in  the  digitized  form  of  the  analytic  signal.  However, 
continuous-time  analytic  signal  generation  is  first  described  for  contextual  purposes. 

The  signal  theory  to  be  relied  upon  is  now  introduced.  Of  particular  importance  is  Euler’s 
identity. 


=cos(a)  +  ysinfa).  (1-1) 

This  provides  a  mathematical  framework  for  interpretation  of  pairs  of  separate  real-valued  signals,  the 
real  and  imaginary  components  of  Equation  (1-1),  as  being  a  single  but  complex-valued  signal.  Here, 
a  represents  some  phase  signal  that  is  a  function  of  either  continuous  or  discrete  time,  i.e.,  either 
a(t)  or  a(nT^ ) ,  and  Ty  is  the  sampling  interval.  Thus  for  an  arbitrary  index  in  time,  n ,  the  in-phase 
signal  is  cos(a)  and  the  quadrature  signal  is  sin(a) .  More  generally,  we  can  multiply  both  sides  of 
Equation  (1-1)  by  a  real-valued  envelope  signal,  |  a  | ,  to  account  for  changes  in  signal  amplitude  over 
time.  As  with  a ,  |  a  |  represents  an  envelope  signal  that  is  a  function  of  either  continuous  or  discrete 
time,  i.e.,  either  |  a(t)  \  or  |  a(nTj, )  | . 
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From  the  modulation  theorem  it  is  known  that  any  real  band-pass  signal,  5(r),  can  be 


represented  as 


s(t)  =1  a(t)  1  •  cos(a(r)) ,  (1-2) 

where  |  a(0  ]  is  the  envelope  of  s(r) ,  and  a(t)  is  the  phase.  If  we  define  the  pre-envelope  of  s(t)  as 


(1-3) 


from  Equation  (1-1)  we  have 

(0  =  s(t)  +  js(t)  =1  a(t)  1  -{cosCorCO)  +  j  sin(a(r))}.  (1-4) 

The  signal  s(t)  =\  a(t)  |  -sinCafr))  is  the  Hilbert  transform  of  s(t)  This  can  be  shown  by  using  the 
mathematical  definition  of  the  Hilbert  transform, 

=  (1-5) 

As  can  be  seen,  this  time-domain  representation  has  little  intuitive  appeal.  The  Hilbert  transform  has  a 
much  more  intuitive  meaning  when  viewed  from  a  frequency-domain  perspective  as  will  be  discussed 
below.  However,  the  time-domain  representation  of  Equation  (1-5)  implies  that  in  order  to  obtain  the 
exact  transform,  the  integral  is  performed  over  all  time.  As  with  an  exact  Founer  transform,  this  is  of 
course,  impractical.  In  practice,  the  transform  is  approximated  with  filters  that  can  actually  be 
implemented.  The  trade-off  areas  in  such  approximations  include  filter  time  delays  and  non-ideal 
frequency  responses  that  are  typical  with  practical  filters. 


*  The  notation  S  is  being  used  here  to  signify  the  Hilbert  transform  of  S  ,  rather  than  the  usual  notation,  S  ;  die  notation  S  will  be  used  to  represent  the 
estimate  of  some  arbitrary  signal,  S  . 
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Figure  1-1.  Representative  band-pass  signal  spectra  from  the  Hilbert  transform  method  of  complex 
envelope  generation;  (a)  the  spectrum  of  the  original  real  signal,  s(t) ,  (b)  the  spectrum  of  j  times  the 
Hilbert  transform,  ?(0,of  s(t),  (c)  the  spectrum  of  the  pre-envelope,  s+{t)  =  s(t)  +  js(t),sa\d  (d) 
the  spectrum  of  the  complex  envelope,  s(t) ,  of  s(t) . 

Figures  1-1  (a)  through  1-1  (d)  help  to  describe  the  Hilbert  transform  process.  In  Figure  1-1 
(a)  is  presented  a  representative  spectrum  of  an  arbitrary  band-pass  signal,  s(t) .  Note  that  the  Fourier 

transform  of  a  real  signal  is  conjugate  symmetric  as  depicted  in  the  figure.  Fourier  transformation  of 
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the  signal  of  Equation  (1-5)  can  be  shown  to  result  in  a  90  degree  phase  shift  in  the  frequency 
components  of  s(t) .  Shown  in  Figure  1-1  (b)  is  the  spectmm  of  the  product  J  -  J(0 .  Multiplication 
of  the  Hilbert  transform  by  the  complex  constant,;,  and  adding  this  to  the  original  signal,  results  in  the 
pre-envelope,  5+  (r) ,  as  described  by  Equation  (1-4). 

Note  that  the  resulting  signal  spectrum  has  no  energy  below  f=0. 
The  utility  of  this  process  is  apparent  through  inspection  of  Equation  (1-4).  With  the  I  signal,  s(t) , 

and  the  Q  signal,  s(t),  one  can  calculate  the  envelope  as  ,  and  the  instantaneous 

phase  as  (r)  =  Arc  tan[?(0/^(0]  •  Alternatively,  the  pre-envelope  can  be  down-converted  m 

frequency  to  0,  by  multiplying  by  ,  where  is  the  center  frequency  of  the  pre-envelope. 

This  results  in  the  complex  envelope. 


sit)  =  s+(,t)’e  (1-6) 

Representing  the  pre-envelope  phase  component  as  a(t)  =  +  ^)(r) + 0 ,  the  phase  of  the  complex 

envelope  is  +  9 ,  where  9  is  some  constant.  The  spectmm  of  the  complex  envelope  is  depicted 
in  Figure  1-1  (d).  hi  practice,  approximations  to  the  complex  envelope  are  more  commonly  used, 
rather  than  pre-envelope  approximations.  As  described  in  Section  3,  however,  this  leads  to 
unnecessary  additional  computations  when  working  with  digitized  signals.  Practical  methods  of 
approximating  the  complex  envelope  can  now  be  presented. 

2.  The  Quadrature  Down-converter 

Also  referred  to  in  the  literature  as  the  quadrature  mixer  or  quadrature  demodulator,  the 
quadrature  down-converter  represents  a'  practical  solution  to  the  problem  of  complex  envelope 
approximation.  The  basic  stracture  of  the  quadrature  down-converter  is  shown  in  Figure  2-1  (a).  It 
consists  of  a  signal  splitter,  and  a  pair  of  mixer  and  low-pass  filter  combinations.  The  intent  of  the 
device  is  to  accomplish  the  elimination  of  the  negative  frequency  components  seen  in  Figure  1-1  (a) 
and  frequency  down-conversion  consistent  with  the  desired  spectmm  of  Figure  1-1  (d).  This  is 
achieved  by  performing  the  down-conversion  first,  then  low-pass  filtering  with  the  appropriate 
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(identical)  low-pass  filter  in  each  channel.  In  the  ideal  case,  if  the  splitters,  mixers,  and  low-pass 
filters  were  ideal  processes  with  perfect  gain  and  phase  responses,  then  a  signal  proportional  to  the 
complex  envelope  is  obtained  by  interpreting  the  device  output  as  /(0  +  (The  proportionality 

constant  factor  of  2  is  easily  accounted  for.)  The  ideal  result  will  be  a  pure  downward  translation  in 
frequency  and  rejection  by  the  low-pass  filters  of  signal  content  that  occurs  in  the  vicinity  of  2/^ . 
This  signal  content  is  expected,  is  due  to  the  mixing  process  itself,  and  is  normally  easily  removed  by 
the  filters.  Hewever,  in  practice,  distortions  arise  in  the  form  of  attenuated  and  phase-delayed  images 
of  the  desired  spectral  components.  In  effect,  the  down-conversion  process  contains  an  undesired  up- 
converted  signal  component  that  is  also  located  at  0  Hz,  creating  the  additive  image  distortion  signal. 


(a) 


Figure  2-1.  The  quadrature  down-converter  method  of  complex  envelope  generation;  (a)  the  physical 
representation,  (b)  the  mathematical  representation. 

The  mathematical  representation  of  the  quadrature  down-conversion  process  is  shown  in 
Figure  2-1  (b).  Comparing  the  physical  and  mathematical  representations,  the  local  oscillator  and 
phase  shifter  processes  are  approximations  to  the  signal  cxp{- .  The  splitter  and  mixer 

processes  attempt  to  achieve  the  same  result  as  complex-valued  multiplication  of  this  locally 
generated  signal  by  the  real  input  signal,  s(t) .  The  pair  of  real-valued  identical  low-pass  filters 
accomplishes  the  rejection  of  out-of-band  frequency  components,  consistent  with  the  rules  of  complex 
arithmetic.  These  rules  imply  the  need  for  a  pair  of  identical  filters,  one  operating  on  the  in-phase 
channel  and  the  other  operating  on  the  quadrature  channel.  In  contrast,  the  mathematical  system 
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requires  a  single  real  low-pass  filter  operating  on  the  complex-valued  signal  output  from  the  complex 
multiplier. 


Figure  2-2.  Representative  band-pass  signal  spectra  from  the  ideal  quadrature  down-converter  method 
of  complex  envelope  generation;  (a)  the  spectrum  of  the  original  real  signal,  5(r) ,  (b)  the  spectrum  of 
the  down-converted  signal,  2sit)  jlTtf.t) ,  (c)  the  spectrum  of  the  complex  envelope,  5(0, 
resulting  from  ideal  low-pass  filtering  and  (d)  the  spectrum  of  the  sampled  complex  envelope, 
5  (nTj) ,  at  the  sampling  interval 
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Once  the  in-phase  and  quadrature  signals  have  been  generated  by  the  quadrature  down- 
converter,  a  pair  of  identical,  synchronized  analog-to-digital  converters  are  used  to  acquire  the  analytic 
signal.  Shown  in  Figure  2-2  are  representative  spectra  for  the  various  steps  in  the  ideal  quadrature 
down-converter,  including  the  sampling  process.  Although  there  are  practical  issues  such  as  gain  and 
phase  imbalance  between  the  channels,  filter  delays  and  local  oscillator  synchronization,  it  can  be  seen 
that  the  quadrature  down-converter  represents  a  viable  method  of  complex  envelope  generation.  In 
fact,  methods  of  characterizing  and  compensating  for  channel  imbalance  have  been  studied  [6-9].  A 
commonly  used  alternative  is  to  first  digitize  the  signal  of  interest  and  then  implement  a  numerical 
(i.e.,  digital)  quadrature  down-converter.  This  leads  to  the  less  known  concept  of  complex  band-pass 
filters. 

3.  Complex  Band-pass  Filters 

The  generalization  of  a  method  of  generating  the  discrete-time  pre-envelope  and  complex 
envelope  signals  from  a  properly  sampled  band-pass  signal  is  now  presented.  Advancements  are 
made  to  the  work  of  Reilly,  Frazer  and  Boashash  [2],  where  the  special  case  which  employs  complex 
FIR  (Finite  Impulse  Response)  band-pass  filtering  has  been  treated.  The  generalization  presented 
herein,  includes  the  extension  to  the  HR  (Infinite  Impulse  Response)  case,  and  the  extension  to  the 
case  where  synchronization  to  a  reference  is  of  concern.  These  extensions  allow  for  further  practical 
application  of  the  complex  band-pass  filtering  process. 

With  the  work  of  [2],  the  concept  of  complex  (band-pass)  filtering  for  generation  of  the 
discrete  analytic  signal  has  been  introduced.  (Related  observations  regarding  the  relationship  of 
symmetric  half-band  low-pass  filters,  Hilbert  transformers,  and  analytic  signal  generation,  have  been 
previously  made  by  Jackson  [10].)  As  the  name  implies,  complex  band-pass  filtering  involves  the 
generation  of  a  pre-envelope  signal,  achieved  via  linear  convolution  of  the  real  input  signal  with  a 
filter  represented  by  an  appropriate  set  of  complex-valued  filter  coefficients.  As  will  be  shown,  these 
coefficients  are  easily  derived  from  the  design  of  a  real  low-pass  filter,  of  appropriate  frequency 
response  characteristics.  This  real  low-pass  filter  will  be  referred  to  herein  as  the  low-pass  filter 
prototype. 

Any  common  method  of  generating  the  coefficients  of  this  low-pass  filter  prototype  may  be 
used,  including  the  popular  Parks-McClellan  algorithm,  and  window-based  designs  [11].  The  same 
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design  considerations  used  in  generating  the  low-pass  filter  of  the  half-band  filter  method  [10]  or  the 
low-pass  filter  pair  of  the  commonly  employed  quadrature  down-conversion  method  [12]  of  complex 
envelope  generation  are  used  in  designing  the  low-pass  prototype.  (For  example,  if  the  cut-off 
frequency  of  the  prototype  filter  is  5/2  Hz,  then  the  bandwidth  of  the  resulting  complex  band-pass 
filter  will  be  5  Hz.)  Once  the  pie-envelope  is  generated  via  complex  band-pass  filtering,  the 
modulation  theorem  indicates  that  one  can  easily  rotate  the  pre-envelope  down  in  frequency  to  a  center 
frequency  of  0,  resulting  in  the  complex  envelope.  Thus  it  can  be  shown  that  the  quadrature  down- 
conversion  and  complex  band-pass  filtering  processes  can  produce  equivalent  results. 

As  in  [2],  this  equivalence  is  leveraged  herein  and  allows  for  the  formulation  of  the  complex 
band-pass  filter  coefficients.  However,  rather  than  being  limited  to  the  FER  case,  the  results  given  in 
Section  3.1  also  apply  to  the  more  general  DR  case.  Interest  in  HR  solutions  is  evidenced  in  [11, 13, 
14].  One  reason  for  this  interest  is  that  it  is  often  the  case  that  FIR  filters  meeting  the  same 
specifications  are  of  substantially  larger  order. 

In  spite  of  the  simplicity  of  the  complex  band-pass  filter,  practical  issues  arise  when 
attempting  to  use  either  this  or  any  equivalent  method  of  analytic  signal  generation.  In  particular, 
situations  often  arise  in  synchronization  applications,  where  filter  time  delays  must  be  taken  into 
account  in  order  to  alleviate  phase  offsets.  In  subsequent  sections  it  is  shown  how  to  avoid  such 
problems.  An  example  of  the  complex  band-pass  filter  method  of  pre-envelope  generation  is  given, 
emphasizing  the  HR  case  and  synchronization  problem. 

3.1  Complex  Band-Pass  Filter  Generation 

Based  on  the  relationship  of  the  quadrature  down-conversion  method  of  complex  envelope 
generation  and  the  complex  band-pass  filter  method  of  pre-envelope  generation,  the  coefficients  of  the 
complex  band-pass  filter  can  be  derived.  Consider  the  low-pass  prototype  filter  impulse  response, 
hpijird ,  and  its  z-transform,  Hip{z).  For  the  ARMA  (Auto-Regressive  Moving-Average)  filter  model, 
which  includes  both  the  FIR  and  HR  cases,  Hip{t)  is  of  the  form 

M 

- •  (3-1) 

k=0 


9 


The  impulse  response,  hip{nTs),  is  deteirnined  by  the  2+M+N  real  coefficients  ftk,  k=0,\,2,...,M  and  r/k. 
k=0,l,2,...J^.  (Here,  is  the  sample  intei-val  in  seconds,  and  the  coiresponding  sample  rate  is 
F,  =  1/  r,  samples  per  second.)  The  filter  response  Htp{z)  corresponds  to  the  response  of  either  of  the 
identical  real  low-pass  filters  used  in  the  quadrature  down-conversion  process. 

The  z-transform  of  a  filter  h{nTs)  is  by  definition 


H{z)=  J,h{nT,)-z-''  . 


(3-2) 


Therefore  when  the  low-pass  filter  prototype  is  appropriately  modulated,  the  resulting  complex 
sequence, 


hip{nTs)e 


+j2n  finT, 


has  the  z-transform 


Thus  for  filters  of  the  form  presented  in  Equation  (3-1)  we  have 


H^iz) 


*=0 _ 


t=0 


fc=0 _ ^ 


i=0 


(3-3) 


(34) 
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Equations  (3-3)  and  (3-4)  imply  that  any  filter  with  the  response  as  indicated  in  Equation  (3-1)  can  be 
circularly  rotated  in  frequency  by  the  amount  fi  Hz,  by  multiplying  the  finite  length  real  coefficient 
sequences  and  by  the  complex  exponential  sequence, 

^+jCl,k  ^^+j2nf,kT,  ^  (3-5) 

The  resulting  complex- valued  coefficients  are  those  of  the  desired  complex  band-pass  filter,  {nT^ ) , 
the  response  of  which  is  given  in  Equation  (3-4).  The  complex  band-pass  filter  numerator  coefficients 
become 


dk  =bk  k  =  0,1,2,... ,M,  (3-6a) 

the  denominator  coefficients  become 

Ck=ak-  ^  =  0,1,2, ...,  .  (3-6b) 

Thus  we  can  transform  a  properly  designed  prototype  low-pass  filter  of  the  quadrature  down- 
converter  with  low-pass  bandwidth  BH,  into  a  complex  band-pass  filter  with  band-pass  bandwidth  B 
and  centered  at  -t-  /  Hz.  Subsequently,  we  can  linearly  convolve  the  input  sequence  with  this  complex 

filter.  The  result  of  this  convolution,  a  positive  pre-envelope,  can  then  be  circularly  rotated  down  in 
frequency  to  yield  the  corresponding  complex  envelope  signal  if  so  desired,  for  e.g.  IM  demodulation 
[15].  This  is  consistent  with  the  results  of  [2],  with  the  important  exception  that  the  method  of 
complex  filter  coefficient  generation  of  Equations  (3-6a)  and  (3-6b)  includes  not  only  the  FIR  case, 
but  also  the  more  general  HR  case.  Note  also  that  we  have  not  limited  ourselves  to  the  case  where  the 
entire  positive  frequency  range  is  of  interest,  as  defined  by  the  pre-envelope  of  the  input  sequence. 
Rather,  the  pre-envelope  of  any  band-pass  constituent  of  the  input  signal  can  be  obtained. 
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3.2  Synchronization  Considerations 

The  rather  simple  looking  results  of  Equations  (3-6a)  and  (3-6b),  are  not  without  associated 
subtleties.  To  better  appreciate  these  subtleties,  we  now  consider  the  application  of  the  complex  band¬ 
pass  filter  method,  to  problems  involving  synchronization.  Without  loss  of  generality,  we  consider  the 
two-component  case  where  a  narrow-band  constituent,  5,  (nT^),  of  the  input  signal, 
sinT^ )  =  (nT, )  +  {nT^ ) ,  is  of  interest.  (The  term  “narrow-band”  is  being  used  relative  to  the 
sample  rate,  /  2 .)  It  is  assumed  that  the  pre-envelope  of  (nT^ )  must  be  extracted,  without  loss 
of  phase  reference  to  s^{nT^).  (A  continuous-time  example  of  such  a  synchronization  problem,  is 
that  of  commercial  FM  stereo  reception.  Subsequent  to  FM  demodulation,  a  pilot  tone  at  19  KHz  is 
used  to  synchronously  demodulate  an  AM-DSB-SC  constituent  at  38  KHz,  which  contains  the  L-R 
audio  component  [16].)  The  approximate  center  frequency  of  is  taken  to  be  £2/  =  27tfiTg 

radians  per  second.  As  it  turns  out,  the  choice  of  the  rotation  sequences  of  Equations  (3-6a)  and  (3- 
6b),  has  an  effect  on  the  ability  to  achieve  synchronization.  In  the  application  of  Equations  (3-6a)  and 
(3-6b),  the  rotation  sequences,  ,  used  in  the  products  shown  were  each  conveniently  assigned  an 
index,  k,  corresponding  to  the  indices  of  the  filter  coefficients. 

This  index  assignment  will  often  suffice,  but  more  generally,  the  indices  must  be  carefully 
selected  to  achieve  a  desired  phase  relationship,  and  need  not  be  integer  valued  in  the  exponent. 
Taking  this  into  account,  we  rewrite  Equations  (3-6a)  and  (3-6b)  as 

d]x-yi  =  ^]x-yi  ■  ,  x  =  y,y  +  l,y  +  2,...,y  +  M  ,  (3-7a) 

x  =  y.y  +  Uy  +  2 . y  +  N.  (3-7b) 

Here,  we  are  using  the  notation  ]jc-j[,  for  the  nearest  integer  value  of  the  quantity  x-y.  In  order  to 
determine  the  proper  values  for  the  index,  y,  and  for  the  offset,  A  ,  we  need  to  consider  the  processing 
involved  in  the  corresponding  quadrature  down-conversion  method  of  pre-envelope  generation  when 
applied  to  the  synchronization  problem.  The  block  diagram  of  such  a  system  for  generating  the  pre¬ 
envelope  of  (nj; ) ,  is  shown  in  Figure  3-1.  This  system  for  generating  the  pre-envelope,  {tiT^ ) , 
of  the  constituent  {nT^ ) ,  has  been  carefully  designed  to  take  into  account  the  effective  time  delay. 
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n^J^,  of  the  low-pass  filter,  hi^(nTJ.  By  pre-determining  and  using  the  sample  delay,  ;z„,  in  the 

system  of  Figure  3-1,  a  pre-envelope  is  generated  that  remains  synchronized  to  the  delayed  input 
signal  constituent,  .  This  further  implies  that  the  phase  relationship  with  the  delayed 

constituent,  s.,([n-n^,]T^),  has  also  been  maintained  and  can  therefore  be  readily  used  to  process 
s.Ain-nJTJ. 


Quadrature  Down-converter 


Figure  3-1.  Pre-envelope  generation  using  the  quadrature  down-converter,  with  phase  relationships  of 
signal  constituents  maintained. 

Analyzing  the  response  of  the  system  of  Figure  3-1,  we  can  determine  the  offset,  A ,  that  will 
yield  the  equivalent  complex  band-pass  filter,  (nT^ ) ,  when  using  Equations  (3-7a)  and  (3-7b).  It 
can  be  shown  that  the  equivalent  complex  band-pass  filter  is  obtained  when  A  =  .  Also,  the 

numerator  filter  coefficients  can  be  made  conjugate  symmetric  by  setting  y  =  -M  /  2 .  This  can  lead  to 
a  lower  number  of  computations  when  implementing  the  filter,  if  the  low-pass  prototype  is  a 
symmetric  FIR  filter. 

As  an  example,  results  from  a  MATLAB-based  simulation  are  presented.  In  this  simulation, 
j, (nT^)  and  s^inT^)  are  both  unity  amplitude  sinusoids  of  frequencies  0.23  and  0.69  Hz 
respectively.  Likewise,  the  phase  angle  of  (jiT^ )  is  three  times  that  of  {nT^ ) .  These  signals  are 
summed  to  form  the  input,  s{nT^ )  =  5,  {nT^ )  +  ^  ^  complex  band-pass  filter  which  is  used  to 
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generate  the  pre-envelope,  s,,{[n-uJT,) ,  of  the  constituent  s.iiiT,).  It  is  then  graphically 
demonstrated  that  the  phase  relationship  with  the  delayed  constituent,  Sjiin-nJT,) ,  has  been 
maintained.  This  is  accomplished  by  tripling  the  phase  and  frequency  of  the  pre-envelope, 
]7-^ )  ^  and  using  the  real  part  of  this  result  as  the  estimate  of  ^2  (t”  ~  ^  •  Samples  of 

both  the  true  constituent,  s,i[n-nJTJ,  and  the  estimate  are  shown  in  Figure  3-2.  (The  starting 

sample  has  been  chosen  well  after  filter  transients  have  subsided.) 

The  prototype  low-pass  filter  is  a  16'*'  order  digital  Butterworth  filter  with  a  specified  cut-off  of 
.09765  Hz  relative  to  F,  =  2  samples  per  second,  resulting  from  the  MATLAB™  “butter.m”  routine 
[17].  In  this  example,  the  pre-envelope  of  a  band-pass  constituent  located  at  /  =0.23  Hz  is  of 
interest.  This  filter  has  been  intentionally  designed  such  that  n,  is  approximately  integer  valued,  and 
in  this  case  n„  =  33  samples.  As  seen  in  Figure  3-2  (a),  when  we  arbitrarily  set  A  =  0,  the  complex 
band-pass  filter  resulting  from  Eqs.  (3-7a)  and  (3-7b)  imparts  an  undesired  phase  on  the  estimate  of 
$2  ([«  -  n„  ]r, ) .  In  contrast  to  this,  when  we  set  A  =  /Iq  =  33 ,  synchronization  between  the  estimated 

and  true  signals  is  achieved  as  shown  in  Figure  3-2  (b). 


Figure  3-2.  (a)  Segments  of  the  estimated  and  true  sequences  of  the  constituent  S2([n  33]r^.), 

with  A  =  0  and  n„  =  33 ;  synchronization  is  not  achieved. 
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Figure  3-2.  (b)  Segments  of  the  estimated  and  true  sequences  of  the  constituent  ^2  ([^  “  33]T^ ) ,  with 
A  =  /!„  =  33  ;  synchronization  is  achieved. 


Figure  3-3.  Example  prototype  low-pass  filter  and  resulting  complex  band-pass  filter  responses;  (a) 
low-pass  filter  magnitude  response,  (b)  low-pass  filter  phase  response, 
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Nomnalized  Angular  Frequency 


Figure  3-3.  (c)  complex  band-pass  filter  magnitude  response,  (d)  compex  band-pass  filter  phase 
response. 


Figure  3-3  shows  the  responses  of  the  prototype  Butterworth  low-pass  filter  and  the  corresponding 
complex  band-pass  filter.  Note  that  the  response  has  been  rotated  in  frequency  to  0.23  Hz. 

33  Example  Application:  Interference  Removal 

The  additional  example  of  maintaining  signal  synchronization  for  the  purpose  of  interference 
cancellation  provides  further  insight  into  potential  applications.  The  cancellation  technique  to  be 
presented  in  this  section  has  been  applied  in  [18]  to  enhance  speaker  identification  in  the  presence  of 
narrow-band  interference.  For  convenience,  a  description  of  the  method  is  also  given  in  this  report. 

Narrow-band  interference  can  be  removed  via  notch  filtering  prior  to  processes  such  as 
demodulation.  For  removal  of  pure  tonal  interference,  the  notch  filter  approach  can  be  very  effective. 
One  can  design  a  notch  filter  with  a  small  bandwidth  such  that  only  a  small  amount  of  the  desired 
signal  is  attenuated  when  the  tone  is  removed.  However,  in  practical  communication  settings, 
interfering  signals  can  have  non-infinitesimal  bandwidths.  A  more  sophisticated  technique  is  required 

to  handle  such  cases. 

The  interference  cancellation  technique/device  proposed  by  the  author  and  presented  herein. 
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will  be  refened  to  as  the  nanow-band  interference  canceller  (NBIC).  The  basic  NBIC  device  is  shown 
in  Figure  3-4.  Although  it  is  somewhat  more  complicated  relative  to  the  notch  filter  method,  it  is  also 
observed  to  be  more  robust  in  scenarios  where  the  interfering  signal  is  angle-modulated.  The  input  to 
the  device. 


XinT,)  =  sinT,)  +  DinT,), 


(3-8) 


Figure  3-4.  A  narrow-band  interference  canceller. 


consists  of  the  sum  of  a  desired  signal  sinT^) ,  and  an  interfering  distortion  signal,  D(nT^) ,  mid  is 
sampled  at  =  l/r,  samples  per  second.  For  the  given  configuration  of  the  NBIC,  it  is  assumed  that 
D(nT^ )  is  much  narrower  in  bandwidth  than  s(nT^) ,  and  can  be  modeled  as 

DinT^ )  =  2|ao  (nT^  )|  •  cos(aro  {nT^ )) ,  (3-9) 

where  ai,{nT,)  =  27tf^nT^  +  4>D(nTJ  +  6o .  Here,  2\a^{nT^^  is  the  slowly  varying  envelope  and  /©  is 
the  center  frequency  of  D(nT^) .  The  angle  modulation  which  is  present  consists  of  <l>D^nT^)  and  the 
constant  (dc)  phase  term,  6d  •  Referring  to  the  figure,  the  function  of  the  complex  band-pass  filter 
(CBPF)  is  to  pass,  with  unity  gain,  frequencies  in  a  narrow  band  about  the  center  frequency  estimate, 
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fo ,  and  to  reject  all  other  frequencies  (both  positive  and  negative).  In  an  initialization  process,  the  set 
of  complex  valued  coefficients  for  the  CBPF  are  generated  based  on  this  center  frequency  estimate. 
To  create  the  CBPF,  Equations  (3-7)  (a)  and  (b)  are  applied,  with  A  =  ,  where  dj  is  the  time  delay 

(in  samples)  imparted  by  the  low-pass  filter  prototype.  The  output  of  the  CBPF  is  an  estimate  of  the 
pre-envelope  of  D(nT^) ,  and  is  represented  as 


Xo(n,T. )  =  1^.  (« J,  )|  •  exp(y  •  {a,  (n.T )  +  e(n,r ) }) ,  (3-10) 

where  n^=n-d^ .  The  envelope  of  the  distortion  estimate  is  \Ao{nT^)\  =  \  0^(117^)  + E{nT^)\.  The 
errors,  E{nT^)  and  e(nTJ  are  an  additive  envelope  error  and  an  additive  phase  error,  respectively. 
These  errors  are  the  combined  result  of  the  component  of  s{nT^)  which  passes  through  the  CBPF, 
and  the  errors  arising  from  the  convolution  of  D{nT,)  with  the  CBPF.  These  effects  are  minimized 
by  proper  selection  of  the  bandwidth  of  the  CBPF.  Nominally,  this  bandwidth  is  set  equal  to  the 
bandwidth  of  D(nT,) .  X^{nT^)  can  be  represented  in  rectangular  notation  as 

XoW  =  X,,inT^)  +  jX,^(nT^) ,  (3-11) 

where 

XoMTJ  =  \Ao(nT^)\-cos{a^(nT^)  +  e(nT^)) ,  (3-12) 

and 

Xo,m  =  \A^(nT,^-sm{a^{nT^)  +  e{nT^)).  (3-13) 

This  representation  facilitates  the  demodulation  processes,  using  the  pair  of  real  signals,  X^inT^) , 
and  Xj^ijtT^) . 
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Before  proceeding  with  the  descriptions  of  the  remaining  processes,  it  is  noted  that  overall,  the 
complex  band-pass  filter,  CBPF,  and  the  low-pass  filters,  LPFl  and  LPF2,  will  impart  a  delay  relative 
to  the  input  signal.  This  fact  is  taken  into  account  in  the  NBIC,  by  designing  these  filters  such  that  a 
total  delay  of  d  samples  is  imparted,  relative  to  the  input.  LPFl  and  LPF2  are  designed  to  impart 
equal  delays  since  the  FM  Demod  and  Integrator  together  impart  a  zero  sample  delay,  as  does  the  AM 
Demod.  Thus  the  delay  of  CBPF  is  d\  samples,  and  the  delays  of  LPFl  and  LPF2  are  both  {d  -d\) 
samples. 

The  AM  demodulation  (AM  Demod)  and  low-pass  filtering  (LPF2)  processes,  are  readily 
accomplished,  and  result  in 


=  \Jxl,([n  - d,]T^)  +  Xl^([n  -  d.W,) 

L  ’  Ju’F2 


=  (3-14) 

The  notation  [•]u.f2  represents  the  low-pass  filtering  process  depicted  in  the  figure  as  LPF2.  From  (3- 
14)  and  the  definition  of  Aj^{nT^) ,  we  have  that 

d,  ([n  -  d]T  )  =  D  ([n  -  d,  ]T  )  +  E({n  -  d,  ]T  )  |]^,, .  (3-15) 

To  be  consistent  with  (3-9),  we  note  the  need  for  a  gain  of  2  as  shown  in  Figure  3-4,  for  restoration  of 
the  envelope.  The  function  of  LPF2  is  to  reject,  to  the  extent  possible,  the  envelope  eiror  component, 

EW  . 

The  FM  demodulation  (FM  Demod)  process  is  implemented  as  a  first-backward-difference 
operation  on  the  recovered  phase  sequence 


=  Arc  tan 


(3-16) 
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This  difference  is  calculated  in  a  modulo-27i  fashion,  often  refeired  to  as  phase  unwrapping.  (More 
conectly,  this  is  refeired  to  as  phase-difference  quadrant  deteiTnination.)  The  Arctan  in  (3-16)  is  the 

four-quadrant  arctangent  of  the  ratio  of  quadrature  and  in-phase  components,  ,  and  ^ oi  •  In  the 
absence  of  phase  cycle-slips,  the  low-pass  filtering  (LPFl)  and  Integiator  processes  result  in  the 
distortion  phase  estimate 

a„([«-dir,)  =  [a„([;!-d,ir,)  +  «(["-^,17:)L,  ■  (3-17) 

The  function  of  LPFl  is  to  reject,  to  the  extent  possible,  the  phase  error  component,  e{nT^) . 
The  Integrator  process,  in  the  case  where  a  backward-difference  FM  demodulator  is  used,  is  simply  an 
accumulator,  with  a  specific  initial  condition.  This  initial  condition  is  set  equal  to  the  initial  phase, 

,  used  in  the  FM  demodulation  process.  (For  simplicity  the  initial  phase  is  set  to  zero.)  Note 
that  for  practical  purposes,  the  integration  process  can  also  be  calculated  modulo-27r,  preventing  large 
accumulation  results.  Thus,  the  function  of  the  Integrator  is  to  invert  the  FM  demodulation  process. 
The  FM  demodulation  process  itself,  has  facilitated  the  use  of  LPFl. 

With  the  envelope  and  phase  estimates  of  D(nTJ  available,  the  cos(-)  and  Multiplier 
processes  are  used  to  generate  the  estimate 

Dan  -  d]T )  =  2a„  ([n  -  d]T )  •  cos{&^  ([«  -  d]T )) .  (3-18) 

Finally,  as  shown,  after  properly  aligning  the  input  with  the  distortion  estimate,  the  distortion  estimate 
is  subtracted  from  the  input  signal,  resulting  in  an  estimate  of  the  desired  signal, 

san-d]Ts)==X  ([;i  -  d]T,  )  -  D([n  -  d]T, ) .  (3-19) 

Thus  the  Time- Alignment  Delay  imparts  a  delay  of  d  samples  on  the  input,  allowing  for  a 
synchronous  subtraction  process. 
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4.  Conclusions 

Generation  of  the  discrete-time  pre-envelope  or  complex  envelope  signals  from  a  sampled 
band-pass  signal,  is  a  fundamental  pre-processing  step  in  modem  signal  processing  systems,  with 
applications  including  sonar,  radar,  audio  and  communication  processing  systems.  The  complex 
band-pass  filter  method  presented  herein  provides  the  system  designer  with  a  simple  means  of  analytic 
signal  generation  which  leverages  existing  real  low-pass  filter  design  tools.  Extensions  are  made  to 
the  cited  previous  work  in  which  the  special  case  of  the  FIR  complex  band-pass  filter  was  presented. 
These  extensions  include  the  more  general  HR  complex  band-pass  filter,  and  applications  where 
synchronization  is  of  concern.  An  example  has  been  presented  which  both  employs  an  DR  complex 
band-pass  filter  and  demonstrates  how  to  maintain  phase  synchronization.  The  concepts  provided  in 
the  example  were  then  built  upon  to  provide  a  more  sophisticated  example  application  of  interference 
cancellation.  The  results  presented  allow  for  further  understanding  and  practical  application  of  the 
complex  band-pass  filtering  process. 

As  presented  in  [4],  the  digital  quadrature  down-converter  represents  an  enhancement  to  the 
traditional  analog  process.  In  particular,  with  appropriately  designed  filters  and  proper  arithmetic 
precision,  the  image  distortion  present  in  analog  implementations  is  kept  within  acceptable  limits  in 
the  digital  implementation.  With  the  work  of  [2],  advantage  is  taken  of  the  equivalence  of  FIR 
complex  band-pass  filters  and  digital  quadrature  down-converters  with  FIR  low-pass  filters.  This 
leads  to  the  elimination  of  the  need  for  complex  sequence  generation  and  multiplication  at  every 
sampling  instant.  Instead,  the  FIR  filters  themselves  are  pre-rotated  in  frequency,  forming  the 
complex  band-pass  filter  used  for  subsequent  convolution.  Designers  of  the  digital  quadrature  down- 
converter  have  cleverly  avoided  complicated  sequence  generation  by  judicious  choice  of  the  sample 
rate  and  center  frequency  of  the  signal  of  interest.  However,  this  limits  the  flexibility  of  the  analytic 
signal  generation  process.  (It  should  be  noted  that  rotation  of  the  center  frequency  of  the  pre-envelope 
signal  to  0  Hz  is  unnecessary;  in  particular,  as  with  envelope  demodulation,  angle  demodulation  does 
not  require  the  complex  envelope  [15].) 

Both  the  FIR  complex  band-pass  filter  process  of  [2]  and  the  digital  quadrature  down- 
converter  of  [4]  address  only  the  case  where  the  entire  positive  frequency  band  is  of  interest.  Li  the 
most  general  case,  any  band-pass  signal  in  the  Nyquist  range  may  be  of  interest.  The  extensions 
presented  in  this  report  allow  for  the  arbitrary  selection  of  center  frequency  and  bandwidth.  More 
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recently,  poly-phase  filter  banks  (see  e.g.,  [19])  are  becoming  popular  for  the  selection  of  band-pass 
signals  that  are  less  than  full  Nyquist  bandwidth.  However,  the  standard  applications  of  such  poly¬ 
phase  methods  are  for  inputs  with  equally  spaced  frequency  channels.  This  and  the  often-used 
additional  constraint  of  perfect  reconstmction,  limits  the  choices  of  filter  responses  available  to  the 
designer.  In  general,  where  no  such  channelization  is  required,  the  complex  band-pass  filter  represents 
an  attractive  alternative  when  signal  fidelity,  bandwidth  and  frequency  selection  flexibility  are  desired. 
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